function y = Gaussian_pulse_delay(t,delay_in)
% Define parameters
A = 1;      % Pulse amplitude (W)
T = 3e-6;      % Pulse period (3 us)
sigma = 0.0000002; % Standard deviation
t0 = T/2;      % Time offset for pulse center
% Calculate current pulse index (n) by floor division
n = floor(t / T);
% Calculate relative time within the current pulse period
t_rel = t - n*T - t0- delay_in;
% Generate Gaussian pulse
y = A * exp(-(t_rel^2) / (2*sigma^2));
end